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Abstract We update our Standard Model predictions for g— 2 of the muon and for the hadronic contributions 
to the running of the QED coupling, Aa^ d (Mf ). Particular emphasis is put on recent changes in the hadronic 
contributions from new data in the 2n channel and from the energy region just below 2 GeV. 
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1 Introduction 

The anomalous magnetic moment of the muon, 
a M = (g — 2)/2, has been the subject of wide interest 
and detailed research. The discrepancy between its 
experimental value as measured by BNL [l| and its 
prediction within the Standard Model (SM) is one of 
the few - if not the only - experimental sign of physics 
beyond the SM (apart from neutrino mixing). This 
has triggered a lot of intense scrutiny of both the ex- 
perimental determination and the theoretical evalua- 
tion of a M . The BNL experiment E821 has achieved 
an impressive precision of 0.5ppm 1], and further im- 
provements may only be reached with the planned 
experiments at Fermilab and J-PARC, see the pre- 
sentations 0, [lj . As discussed below, the SM predic- 
tion relies heavily on the experimental information of 
the measured hadronic cross sections at low energies. 
During the last 15 years, the SM prediction has gone 
through several phases of improvement and consol- 
idation and has now, for the first time, reached an 
accuracy even slightly better than the experimental 
value. This is mainly due to the big efforts to measure 
the hadronic cross sections with increasing accuracy 
and the progress of various groups working on the 
data-driven evaluation of the hadronic contributions, 
which are in fairly good agreement (though further 
improvements are foreseen as will be discussed briefly 
in Section |4j . In this article we will concentrate on 
the main changes in the hadronic contributions to a M , 
updating the works 0, |B| ■ 

Similarly to g — 2, the theoretical uncertainties 
in the running of the QED coupling, a(q 2 ) 1 are 



completely dominated by the hadronic contributions, 
Aa^ d (q 2 ). Recall a[M%) is the least well known of 
the set of precision observables [G F ,M z ,a(M^)), so 
its error is a limiting factor in the electroweak fits of 
the SM as performed e.g. by the LEP Electroweak 
Working Group. 

2 Standard Model prediction of g — 2 

The anomalous magnetic moment of the muon re- 
ceives contributions from all sectors of the SM. The 
QED and electroweak (EW) corrections can be cal- 
culated within perturbation theory and are, through 
many impressive works, well under control. In the 
compilation for a^ M presented here we use a^ ED = 
116584718.08(15)10- 11 [IB and a% w = (154±2)1CT 11 
8], as e.g. reviewed in |J The hadronic contribu- 
tions can not be reliably calculated in perturbative 
QCD (pQCD) as the loop-integrals are dominated by 
low momentum transfer, i.e. the non-perturbative 
region of QCD. They are typically divided into the 
leading (LO) and higher-order (HO) vacuum polari- 
sation (VP), and the so-called light-by-light contribu- 
tions, which are also subleading: a^ ad = a h ad . LOVP + 
had,HOVP + a had,i-by-i < while the V p induced cor- 

/J, f-L 

rections can be calculated with methods based on 
dispersion relations and using experimentally mea- 
sured hadronic cross sections as input, the light-by- 
light scattering contributions can be estimated only 
using modelso The results from different groups vary 
considerably, both w.r.t. the mean value and the er- 
ror. For a review presented at this conference see 10 1. 
In the SM prediction of g—2 presented here we use the 
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value a^i-by- 1 = (10.5 ±2.6) • KT 10 , which has been 
obtained in [ll[ as a combination of results based on 
different modelsjj] In the following we will discuss in 
more detail recent changes in the VP contributions. 

2.1 Hadronic VP contributions 

The LO hadronic VP contributions are calculated 
using the dispersion integral 

1 



had.LOVP _ 



4vT 3 



(1) 



where K(s) = • (0.4. . . 1) is a known kernel func- 
tion giving highest weight to lowest energies ^fs, and 
f7 had( s ) i s the hadronic cross section for e + e~ — > 7* — > 
hadrons (+7) . The superscript indicates that the 
'undressed' cross section must be used, i.e. the cross 
section without VP effects in the virtual photon, but 
including final state radiation (FSR) of photons. To 
arrive at the best compilation for c had , at low ener- 
gies (y/s < 2 GeV) about 24 hadronic channels (ex- 
clusive final states) have to be summed, and in each 
channel the data from different experiments have to 
be combined. At intermediate energies tr h ad is mea- 
sured inclusively. Perturbative QCD can only be used 
away from resonances, and (most) data driven anal- 
yses use pQCD only for energies above the open bb 
threshold. For details of the data input, their treat- 
ment w.r.t. radiative corrections and the data com- 
bination through a non-linear Xmin fit, as used in the 
work reported here, see 0, B|- Note that there are 
uncertainties w.r.t. the correct application of radia- 
tive corrections (undressing of VP, possible addition 
of neglected photon FSR) to the data, especially in 
the case of older data sets. In our analysis this leads 
to the assignment of a separate error due to radiative 
corrections, Sa h ^' vp+FSR ~ 1.8 • lO" 10 , which alone is 
about ten times bigger than the uncertainty of the 
electroweak contributions a^ w . 

2.1.1 Recent changes in the 2tt channel 

More than 70% of a^ ad is coming from the p 2tt 
channel. Following older experiments, the cross sec- 
tion e + e~ — > tt + tt~ has been measured in recent years 
with increasing accuracy by the Novosibirsk experi- 
ments CMD-2 and SND, see e.g. [l4] for a short re- 
view of their results. Figure [1] displays the impressive 
agreement of the most recent data sets from CMD-2 
and SND, together with recent data from KLOE 15] 
obtained with the method of Radiative Return (see 
for a detailed review of the method and its ap- 
plication, and [17[ for the very new KLOE analysis 



of the 2ir channel presented at the PhiPsi09 confer- 
ence but not yet available for public use). The grey 
band shows the result of our data combination in this 
channel including also older data, but excluding the 
KLOE data, see the discussion below. 



Fit (w/o KLOEI 
Mean of fit (w/o KLOE) 
Mean of fit (all sets 

KLOE (08) t 
CMD-2 (06) t 
SND (04. re-anal.) ^ 
CMD-2 (02. re-anal.) t 




Energy [GeV] 

Fig. 1. Most important data in the 2ir channel 
and fits as described in the text. 
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Fig. 2. Low energy region close to the 2n threshold. 

The lowest energy region close to the 27r thresh- 
old is displayed in Fig. [2] The recent CMD-2 data 
represented by (red) triangles and error bars clearly 
demonstrate the improvement from this single data 
set alone, with the fit in this region being dominated 
by these data. Figure [3] shows an enlargement of the 
p-uj interference region which is now very well mapped 
out by the consistent data sets. However, Fig. [3] also 
shows that the KLOE data are undershooting the 
combination of the other data in this region, while 
typically being higher than other data at lower ener- 
gies as can be seen by the (green) solid line in Fig. [TJ 
This apparent difference in shape is highlighted in 
Fig. |4j which shows the normalised difference of the 
KLOE cross section and the combination of the other 



T Note that this is slightly different from (though certainly compatible with) the value Ou a ' 
m and discussed in lid. 
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: (116±40)-10 -11 as obtained 
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2tt data by the square markers, while the band dis- 
plays the size of the error of the compilation with- 
out KLOE. (The bands displayed in the figures are 
obtained from the diagonal elements of the fit's full 
covariance matrix.) 




0.77 0.78 0.79 

Energy [GeV] 



Fig. 3. p-uj interference region in the 2ir channel. 
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Fig. 4. Normalised difference of the KLOE 
data [lfj and the data compilation excluding 
KLOE; the band represents the error of the 
compilation. 

One should keep in mind that the KLOE data 
are obtained via the method of Radiative Return 
at fixed collider centre-of-mass energy, whereas the 
other data are measured via the traditional method 
of energy scan by adjusting the e + e~ beam ener- 
gies. Hence Monte Carlo simulation tools including 
radiative corrections and also the systematic effects 
are completely different between the two approaches. 
Up to now it is not clear what causes the differ- 
ence in the shapes of the 2tt data, and more stud- 
ies are underway to clarify the situation. Unfortu- 
nately this difference in shape prevents us from in- 
cluding the KLOE data in a straightforward way in 



the non-linear x^ in Note that if we calculate 

the contribution to g — 2 from the KLOE data alone, 



we obtain 



, KLOE _ 



(384.16±3.47)-f0- 10 , in per- 



fect agreement with the result of the integral over 
our compilation of all 2tt data without KLOE (but 
in the range of the KLOE data), for which we get 
^w/out KLOE = (384.12 ±2.51) ■ lCT 10 . We therefore 
use the procedure adopted already in 0] where earlier 
KLOE data (now superseded by [15j) were included 
after integration, by calculating a weighted average 
with the integral over the data compilation without 
KLOE, and not performing a point-by-point combi- 
nation. 

BaBar has also published their first measurement 
of the 27T channel based on Radiative Return [3| (see 
also 19]), finding some discrepancies with KLOE. We 
have not included the new BaBar data in the analysis 
presented here as they were not available for public 



ilEl for 



use at the time of the conference, but see 
an analysis which includes them. 
2.1.2 Energy region below 2 GeV 

Important changes in the data input have also 
happened in the region between 1.4 and 2 GeV. This 
region is particularly difficult, as a growing number 
of multi-hadron exclusive final states becomes ac- 
cessible and has to be included to obtain cr ha d with 
good accuracy. However, these energies are above 
the reach of the Novosibirsk machine^, and the qual- 
ity of the available data was not very good. Alter- 
natively, one can rely on inclusive R measurements, 
but for this only rather old and not very precise 
data are available. This situation has changed with 
BaBar measuring, through Radiative Return, many 
channels with higher accuracy than earlier experi- 
ments. These include new data for 2/T + 27r~ |23j |. 
K°irK F 



2tt+2tt-tt°, K+K-ir+w-w 



K+K- 

2tt + 2tt-t] 0, 27r + 27r-27r° [26| used for our updated 
analysis. Figures [S] - [7] exemplify the influence of the 
new BaBar data. The new data are not always in 
good agreement with other sets; in such cases the fit 
has a poor quality and we scale up the error of the 
channel's contribution by yj Xmm/d-O.f. (e.g. in the 
27T+27r-27r channel xL n / d - oi - = 2 - 7 -) 
2.1.3 Sum rule analysis 

Note that in previous g—2 analyses 0,01 we found 
a discrepancy between using the available inclusive 
data in the region 1.4 to 2 GeV and adding the ex- 
clusive channels. The inclusive data were lower than 



+ The fit allows a readjustment of the overall normalisation of the data sets within their systematic errors. Including the KLOE 
data would lead to a bad Xmin/d-°-f- an< ^ unnatural normalisation effects pulling the fit upward, see § for a detailed discussion. 
3 This will change with the current upgrade, see lii. 
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the sum of the exclusive channels, though they were 
found to be similar in shape. 



* 1 
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Fig. 5. KgK-K channel with improvement due 
to recent BaBar data [24|. 



over the exclusive for energies y/s = 1.43... 2 GeV. 
Since then the situation has changed: the hadronic 
data has changed slightly, being lower at low energies 
and also at energies above 2 GeV. Also, with the in- 
clusion of the recent BaBar data the sum of the exclu- 
sive channels in the region 1.4 to 2 GeV has become 
slightly lower and more accurate than before. Our up- 
dated sum rule analysis is summarised in Fig. [3] dif- 
ferent sum rules based on pQCD are made to match 
the corresponding sum rule integrals over the data by 
fitting for a a as a free parameter (see [1] for details). 
It is clear that the new sum over exclusive channels 
is more accurate than the old inclusive data and also 
more compatible with the predictions based on pQCD 
with a world average value of a a . We therefore are 
now combining the results from the inclusive and the 
sum over exclusive data in this energy region0 




Fig. 6. 2tv + 2tt channel. 




Fig. 7. 27T + 27r~27r° channel. 

In Q we performed a QCD sum rule analysis, con- 
cluding that the inclusive data are more compatible 
with pQCD and the world average of ot s , and there- 
fore chose to use the inclusive data instead of the sum 
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Fig. 8. Results for different sum-rules [H] trans- 
lated into a prediction of a s ; the band shows 
the world average of a. (Ml). 

2.1.4 Other changes and result for a^ d,VP 

In addition to the important changes discussed 
above, compared to [4|, we have also included new 
data in other channels: K + K~ from CMD-2 [i^] and 
SND Q, K° S K° L from SND @, vt+tt-tt from CMD- 
2 0, W7T° from KLOE EH, and inclusive R data at 

and 



higher energies above 2 GeV from BES [32J, [3 
CLEO [34]. Figure H shows the recent BES data to- 
gether with the fit of all inclusive data in this region 
and the prediction from pQCD. While the contribu- 
tion to g— 2 is significantly smaller than previously, it 
is obvious that pQCD, which is in perfect agreement 



^Note that the sum over exclusive channels stills requires, due to the lack of experimental information, the use of iso-spin 
relations for unknown channels, which in turn results in a large error from the poorly known KKtttt channel. 
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with the three latest BESII (09) data points |33[, is 
still somewhat lower than the fit for energies below 3 
GeV. 




Fig. 9. Recent data from BES together with the 
fit of all inclusive data compared to pQCD in 
the energy region above 2 GeV. 

Numerically the changes to a h ad > LOVP from the 
different energy regions amount to (units of 10~ 10 , 
compared to results from fl): -0.76 (0.32-1.43 GeV, 
low energy exclusive channels), +2.10 (1.43 — 2 GeV, 
results from inclusive and sum over exclusive data 
combined), —1.35 (2 — 11.09 GeV, higher energy in- 
clusive data). This accidentally leads to a near per- 
fect cancellation of the shifts, with the total result 
a had,Lovp = (689.4±3.6 cxp ±1.8 rad )-10- 10 . The first 
error is coming from the statistical and systematic 
error of the data, whereas the second error is our 
estimate of the uncertainty in the radiative correc- 
tions as mentioned above. However, compared to 
our earlier result from 4], aJ; ad ' LOVP (HMNT06) = 
(689.4±4.2 cxp ±1.8r ad )-10- 10 , there is a further reduc- 
tion in the error. 

The higher-order VP induced contributions can 
also be calculated using a dispersion integral, see jjjj 
for details. Our updated result is nearly unchanged 
and reads a |; ad > HOVP = (-9.79±0.06 cxp ±0.03 rad >10- 10 . 

2.2 SM predictions compared to the BNL 
measurement 

Combining the QED, EW and hadronic contribu- 
tions as discussed above, we arrive at our SM predic- 
tion of the anomalous magnetic moment of the muon, 



'(HLMNT09) = (116 591 773±48)-10" 



(2) 



Due to a small shift of Codata's muon to proton mag- 
netic ratio A the experimental value for a M from BNL 
is now a^ xp = 116 592 089(63)-10- n 0. This results 



corresponds to a 4 a discrepancy. The situation is 
displayed in Fig. 1101 where we also show our previous 
result and the most recent results from Jegerlehner 
and Nyffeler 1131 and three different predictions from 
Davier et al. [20|, |2l|, |35 1 : While their e + e~ based re- 



sult without the new BaBar 2ir (labelled 'w/o BaBar 
(09)' in the figure) data agrees very well with our 
evaluation, the result including these new data lead 
to a shift upwards, but still compatible with the other 
e + e~ based result. The third result employs, in addi- 
tion to e + e~ data, also the use of r spectral function 
data from ALEPH, OPAL, CLEO and BELLE. 



HMNT (06) 
JN (09) 

Davier et al, x (09) 



w/o BaBar (09) 
w/ BaBar (09; 



Davier et al, e + e" 
HLMNT (09) 
— experiment — 



BNL 

BNL (new from shift in A.) 

iiiiIiiiiIiiiiIiiiiIiii 

a^ SM x 10 10 - 11659000 



i I i i i i I i i i i I i i i i I i i i ill i 

170 180 190 200 210 



Fig. 10. Comparison of recent predictions for 
g — 2 compared to the BNL measurement. 

To translate the charged current induced hadronic 
r decay data into the required spectral functions 
requires the application of isospin breaking correc- 
tions, which can only be predicted in a model de- 
pendent way. Earlier r based results from Davier et 
al. were incompatible with e + e~ based results, but 
with a re-evaluation of the isospin breaking correc- 
tions 



21 



35 



they find the result displayed in 
Fig. [TU] (labelled V (09)') which is now marginally 
consistent. These findings were discussed controver- 
sially at the PhiPsi09 conference; Benayoun presented 
an alternative approach based on Hidden Local Sym- 
metry and dynamical (p,aj,r/>) mixing. With this and 
with a global fit he gets consistency of the r spec- 
tral function with the e + e~ data and an improved 
2tt contribution to g — 2 [37] • Keeping in mind that 
other studies obtained a very large uncertainty of the 
isospin breaking corrections (see e.g. [38]) we believe 
that, until this issue is better understood, the pre- 
dictions of g — 2 based on e + e~ data alone are more 



in a difference a EXP 



(31.6±7.9)-10- 10 which reliable. 
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3 Running QED coupling a(q 2 ) 

Based on the same data compilation, we also pre- 
dict the hadronic contributions to the running of the 
QED coupling, a(q 2 ) = a/(l — Aa ]cp — Aa had ). The 
hadronic VP is important for many studies where 
high accuracy radiative corrections are required, as 
is the case for the hadronic contributions to g— 2, see 
16j for a very recent review. Of particular impor- 
tance as input in EW precision fits is the quantity 
Aa[^' d (M|), the hadronic contributions to the run- 
ning a from five-flavours (the contribution from the 
top quark is usually added using pQCD). Our up- 
dated value is Ao^Mf) = 0.02760 ±0.00015. This 
is slightly higher and significantly more accurate than 
the number from Burkhardt and Pietrzyk used as de- 
fault by the LEP EW Working Group, and leads e.g. 
to a lower preferred Higgs mass and a lower upper 
limit from the famous 'Blue Band Plot'. 



4 Summary and outlook 

We have given an update of the SM prediction 
of g — 2, emphasising recent developments for the 
hadronic contributions. With the e + e~ based anal- 
ysis presented here we obtain a 4 a discrepancy be- 
tween the experimentally measured value of a M and 
its SM prediction. Slightly less significant discrepan- 



cies are reported by other groups, depending on the 
data used and the details of the analyses. However, 
for all e + e~ based analyses the discrepancy is about 
3 — 4 a and standing all scrutiny. 

For the future, further improvements of the SM 
prediction will be possible. The method of Radiative 
Return has proven extremely powerful and will be 
leading to many more results. In addition to the al- 
ready reported new 27r data from KLOE 17J , further 
analyses are ongoing, and there are exciting prospects 
with KLOE2 [39(. There is also a rich programme go- 
ing on at BELLE, and the possibility of SuperBELLE 
at the horizon. CMD-3 and SND at VEPP-2000 
currently commissioned in Novosibirsk are aiming at 
largely improving the exclusive measurements in the 
region below 2 GeV, whereas BESIII at BEPCII will 
cover higher energies. With all these developments 
the error of the SM prediction of g — 2 is expected to 
shrink even further, so that in the future the light-by- 
light contributions may eventually become the limit- 
ing factor. Obviously there is an extremely strong 
case for new the g — 2 experiments as planned by the 
2 Collaboration for Fermilab |2J and at J-PARC 
With this, g — 2 will become even more powerful 
in establishing and constraining physics beyond the 
Standard Model. 
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